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ABSTRACT 

We outline the possible physical processes, associated timescales, and energetics 
that could lead to the production of pulsars, jets, asymmetric supernovae, and weak 
7-ray bursts in routine circumstances and to a 10 16 G magnetar and perhaps stronger 
7-ray burst in more extreme circumstances in the collapse of the bare core of a massive 
star. The production of a LeBlanc- Wilson MHD jet could provide an asymmetric 
supernova and result in a weak 7-ray burst when the jet accelerates down the stellar 
density gradient of a hydrogen-poor photosphere. The matter-dominated jet would be 
formed promptly, but requires 5 to 10 s to reach the surface of the progenitor of a Type 
Ib/c supernova. During this time, the newly-born neutron star could contract, spin 
up, and wind up field lines or turn on an a — dynamo. In addition, the light cylinder 
will contract from a radius large compared to the Alfven radius to a size comparable 
to that of the neutron star. This will disrupt the structure of any organized dipole 
field and promote the generation of ultrarelativistic MHD Waves (UMHDW) at high 
density and Large Amplitude Electromagnetic Waves (LAEMW) at low density. The 
generation of these waves would be delayed by the cooling time of the neutron star 
~ 5 to 10 seconds, but the propagation time is short so the UMHDW could arrive at 
the surface at about the same time as the matter jet. In the density gradient of the 
star and the matter jet, the intense flux of UMHDW and LAEMW could drive shocks, 
generate pions by proton-proton collision, or create electron/positron pairs depending 
on the circumstances. The UMHDW and LAEMW could influence the dynamics of 
the explosion and might also tend to flow out the rotation axis to produce a collimated 
7-ray burst. 



Subject headings: supernovae: general — gamma rays: bursts — pulsars:general — ISM: 
jets and outflows 
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1. Introduction 

Recent evidence has given support for the idea that the core-collapse process is intrinsically 
strongly asymmetric. The spectra of Type II and Type Ib/c supernovae are significantly polarized 
indicating asymmetric envelopes (Mendez et al. 1988; Jeffrey 1991; Hoflich 1991; Trammel et 
al. 1993; Wang et al. 1996; Tran et al. 1997; Leonard et al. 1999). The degree of polarization 
tends to vary inversely with the mass of the hydrogen envelope (Wang, Wheeler & Hoflich 2000). 
Pulsars are observed with high velocities, up to 1000 km s _1 (Strom et al. 1995). Observations of 
SN 1987A showed that radioactive material rapidly mixed out to the hydrogen-rich layers (Lucy 
1988; Sunyaev et al. 1987; Tueller et al. 1991). Cas A shows rapidly moving oxygen-rich matter 
outside the nominal boundary of the remnant (Fesen & Gunderson, 1996) and evidence for two 
oppositely directed jets of high- velocity material (Reed, Hester, & Winkler 1999). High velocity 
"bullets" of matter have been observed in the Vela supernova remnant (Taylor et al. 1993). 
Other evidence shows that soft gamma-ray repeaters arise in very strongly magnetized neutron 
stars, "magnetars," with dipole fields in the range 10 15 G (Kouveliotou et al. 1998; Thompson 
et al. 1999). Theoretical models have shown that matter-dominated jets can cause supernova 
explosions (Khokhlov et al. 1999) and that a strongly asymmetric explosion can account for the 
outward mixing of 56 Ni in SN 1987A (Nagataki 1999). In addition, much attention has recently 
been paid to the issue of collimation of 7-ray bursts and the associated affect on energetics and 
observable properties in the context of high-redshift events (Kulkarni et al. 1999; Rhoads, 1999; 
Sari, Piran & Halpern 1999; Lamb 1999; Chevalier & Li 1999a,b) and the proposed correspondence 
of SN 1998bw with a much lower energy event (Galama et al. 1998). In this paper, we explore the 
physics that might unite all these areas. 

The discovery of optical afterglows of 7-ray bursts has raised the estimates of the maximum 
isotropic 7-ray emission to unprecedented values, ~ 3 x 10 54 ergs for GRB 990123 (Kulkarni, et 
al. 1999). This has brought a new focus on the liklihood that the prompt 7-ray bursts and the 
afterglows are collimated to various extents as well as Doppler boosted and "beamed" (Rhoades 
1997, 1999; Sari, Piran k Halpern 1999; Kulkarni et al. 1999; Stanek et al. 1999; Harrison et 
al. 1999). Although the evidence is still preliminary, collimation factors of AJ7/4-7T ^ 0.01 have 
been derived from the decline in some afterglow light curves. Wang & Wheeler (1998), Nakamura 
(1998) and Cen (1998) have pointed out that if the collimation were strong enough, the energetics 
might be provided by supernova-like energies. SN 1998bw and GRB 980425 provided a different 
perspective by suggesting that some 7-ray bursts are directly associated with some supernovae 
(Galama et al. 1998; see also Bloom et al. 1999; Reichart 1999; Germany et al. 1999; Galama et 
al. 1999; Wheeler 1999). If this association is correct, the "isotropic" 7-ray energy of GRB 980425 
is only ~ 10 48 ergs. Some models of SN 1998bw invoked especially large kinetic energies, in 
excess of 10 52 ergs in spherically-symmetric models, to account for the bright light curve and 
high velocities (Iwamoto et al. 1998; Woosely, Eastman & Schmidt 1998), and others took note 
of the measured polarization and chemical structure to suggest that strongly asymmetric models 
could account for the observations with more "normal" energies (Hoflich, Wheeler & Wang 1999; 



-3- 



Danziger et al. 1999). It is not clear that either class of models can naturally produce a 7-ray 
burst of ~ 10 48 ergs. The former models have been referred to as "hypernova" models, although 
that term was originally introduced (Paczyhski 1998) to mean the generic high energy events 
associated with the afterglows of the classical 7-ray bursts at large distance. 

Many of the popular models for creating 7-ray bursts are based on binary neutron stars 
(Paczyhski 1986) or accretion onto black holes (Paczyhski 1991; Woosley 1993). The latter, in 
particular, are popular because there is no limit, in principle, to the mass of the black hole and 
hence, again in principle, the energy that can be extracted. All the models struggle with the 
mechanism for turning the large energies into 7-rays. There is no reason why the energy flux from 
a black hole source cannot be collimated, as witness the jets from active galactic nuclei and some 
binary black hole sources. This could reduce the energy requirements even in black hole models. 
Some models explicitly include this collimation (MacFadyen & Woosely 1999). We note that even 
the superluminal jets from AGN and blazars have a maximum bulk Lorentz factor of about 10, 
while the 7-ray bursts require T <; 100 (Krolik & Pier 1991; Baring & Harding 1997: Piran 1999, 
and references therein). This may suggest that a qualitatively different mechanism is needed to 
generate the cosmic 7-ray bursts despite the attractive possibilities for black hole models. 

Here we will attempt to see how far a more conservative model can go, both to produce 
asymmetric supernovae and perhaps to generate 7-ray bursts, by considering the effects of a 
newly-born pulsar. This is not a new idea (Ostriker & Gunn 1971; Bisnovatyi-Kogan 1971; 
Bisnovatyi-Kogan et al. 1975), but it is worth reconsidering in the context of the polarization 
of core-collapse supernovae, the growing evidence for "magnetars," the growing understanding 
of 7-ray bursts and their afterglows, the likely association of SN 1998bw with GRB 980425, and 
the interesting possibility that supernovae can be induced by energetic jets arising from their 
cores. All these processes demand or suggest strong asymmetries. Questions arise as to whether 
these issues are related, whether 7-ray bursts of observed properties could be produced in this 
context, whether there are more than one type of 7-ray burst mechanism, and whether supernovae 
may, in some circumstances, produce the 7-ray bursts seen at high redshift. It is very plausible 
that the formation of a neutron star will engender collimated flow, including Large Amplitude 
Electromagnetic Waves (LAEMW: Usov, 1992, 1994; Blackman & Yi 1998). Here we will show 
that it is plausible that the properties of SN 1998bw, including a weak 7-ray burst, may be 
generated by the acceleration of a collimated shock down a density gradient and possible that 
a 7-ray burst visible at cosmic distances could be produced under some circumstances. Certain 
aspects of our model are similar to those of Usov (1992, 1994), but he considered accretion-induced 
collapse and we explicitly investigate the context of a core collapse in the core of a massive star. 
Some of the ideas we explore here are also presented by Nakamura (1998). 

The purpose of this paper is to outline the basic time scales, energetics, and relevant physical 
processes in order to define areas that need more quantitative work. Key length and time scales 
and a summary of the core collapse ambiance are given in §2 and §3, respectively. The physics 
of the proto-neutron star phase including the generation of the magnetic field and associated 
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torquing of the surrounding plasma is presented in §4. The effect of an axial jet associated with 
the collapse process is outlined in §5. The important phase when the neutron star cools, contracts, 
and spins more rapidly is discussed in §6, and the manner in which the energy associated with 
matter and radiation jets could be propagated outward is given in §7. Discussion and conclusions 
are given in §8. 



2. Basic Length and Time Scales 

In the discussion to follow there are a number of key length and time scales. Among these 
are the radius of the star, i? s tar — 10 8 km for a red supergiant and Ru c — 2 x 10 5 km for a helium 
core. The inner iron core that collapses to form a neutron star has a typical radius, Rp c ~ 4 x 10 3 
km. The dynamical or sound crossing time is 

R 3/2 

rd ^ 20S (M/M Q )V 2 ' « 

where R§ is the radius in units of 10 5 km, and the light travel time is 

R 

riight = -^0.3si? 5 . (2) 

The dynamical time to emerge from a red giant or even a light crossing time, from 10 3 to 10 5 s, is 
not likely to be associated with observed 7-ray bursts. For a stripped helium core, the progenitor 
of a Type lb or Type Ic supernova, a range of timescales from 1 to 10 s is possible, in the absence 
of Doppler boosting. Note that the breakout of a collimated jet might involve an area considerably 
less than the surface area of the helium core and an appropriately smaller time scale. Other 
relevant length scales are the radius of the light cylinder, 

Rlc = £ = f- ^ 50 km (—) , (3) 
iZ 2tt \ms/ 

where fi is the rotational frequency and P the rotational period of the neutron star, and the 
Alfven radius at which magnetic pressure is balanced by the ram pressure, e.g., 



k****- (4) 



which, with B ~ B NS (jjfj) , is 



B \V3/ p x-l/6/ y v-1/3 



During the thermal contraction of the proto- neutron star, P, Rlc, and Ra will all change 
significantly. For illustration, we assume that the rotation period decreases from that of the 
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proto-neutron star, Ppns — 25 ms, to that of the neutron star, P^s — 1 ms. As a result, the light 
cylinder will contract from Rlc — 10 3 km, beyond the radius of the stalled shock and comparable 
in size to the original iron core, to Rlc — 50 km, comparable to the radius of the neutron star 
and well within the stalled shock. The light cylinder will contract from beyond Ra to significantly 
less than Ra- This has critical implications for angular momentum transport and the generation 
of radiation. 



3. Context of Core Collapse 

We know that pulsars arise from core-collapse. We know that some core collapse events are 
directly associated with supernovae, e.g., the Crab Nebula, SN 1987A. Circumstantial evidence 
suggests that supernovae of Type II and Ib/c result from core collapse based on correlations with 
spiral arms and H II regions in spiral galaxies and on the nebular spectra that are consistent 
with the cores of massive stars. Polarization data suggests that all these events are significantly 
asymmetric and that the asymmetry is stronger for explosions with smaller hydrogen envelopes 
(Wang, Wheeler & Hoflich 1997; Wang & Wheeler 1998; Hoflich, Wheeler & Wang 1999; Leonard 
et al. 1999; Wang, Wheeler & Hoflich 2000). Explosions induced by jets from within the inner core 
can plausibly account for this asymmetry (Khokhlov et al. 1999). The imprint of such jets and 
the possibility to make 7-ray bursts is largest in massive stars from which the hydrogen envelope 
has been stripped by winds or binary mass transfer. In the following we will thus concentrate on 
core collapse explosions and neutron star formation in hydrogen and helium-deficient Type lb and 
Type Ic supernovae. Such supernova progenitors must be surrounded by substantial mass from 
the previous mass loss stages, and this will affect the production and evolution of jets and 7-ray 
bursts. 

Here we will adopt a generic picture of core collapse in an intermediate mass star of main 
sequence mass ~ 15 - 25 M with length scales as given in §2. The iron core undergoes instability 
and collapses in a time of ~ 1 s. The collapse involves a homologous collapse of the iron core in 
which the density increases monotonically inward. The outer parts of the stellar core composed 
of silicon, oxygen, carbon, etc, have lower densities and longer free-fall times. Over the times of 
interest here, ~ 10 s, these outer layers will "hover" as the collapsing iron core succeeds or fails to 
generate an explosion. 

After core bounce, a proto-neutron star (PNS) forms with a radius of Rpns — 50 km and a 
shock is formed that stalls at a radius of R s h — 200 km. Over a cooling time t coo \ ~ 5 to 10 s, 
the proto-neutron star de-leptonizes by neutrino emission, cools, and contracts to form the final 
neutron star structure (Burrows & Lattimer 1986). If the neutron star is rotating and magnetized, 
this cooling phase will be associated with contraction, spin-up, and amplification of the magnetic 
field. These changes can significantly alter the physics associated with the neutron star and its 
interaction with its surroundings and hence the explosion process itself. 
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4. The Proto-Neutron Star Phase 

Right after collapse, the hot proto- neutron star undergoes neutrino-driven convection with the 
characteristic convective overturn time scale T conv — 1 ms F 39 1/3 where F 39 = F/10 39 erg s 1 is the 
neutrino flux at the base of the convection (Duncan & Thompson 1992). The hot proto-neutron 
star has a radius ~ 50 km that gradually shrinks on the cooling time scale ~ 5 to 10 s. This time 
scale is coincidently, but significantly, similar to the sound crossing time of the helium core, and 
hence about the time required for a supersonic (but not relativistic) jet to penetrate the outer 
mantle. There are thus two potential mechanisms to create a 7-ray burst about 10 s after collapse. 
A 7-ray burst could be generated as the jet accelerates down the density gradient at the boundary 
of the core. Alternatively, a strong flow of Poynting flux might be generated on the neutron-star 
cooling time as the neutron star spins up. This might lead to an alternate mechanism of 7-ray 
burst formation on a similar timescale (Nakamura 1998). These mechanisms for generating 7-ray 
bursts will be discussed in §7. 

It is hard to estimate the spin period of the proto-neutron star when it first forms. We 
assume that any 7-ray burst stimulated by Poynting from the neutron star will occur near the 
end of the cooling, contraction stage when the neutron star spins with a period Pns — 1 ms - As 
we will see below, a spin period approaching Pns — 1 ms is necessary (if not sufficient) to power 
a classical 7-ray burst. A slower spin or shock breakout associated with a jet might generate the 
type of 7-ray burst observed in SN 1998bw/GRB 980425 as we will also show below. For purposes 
of illustration, we then adopt a spin period of the proto-neutron star to be that which will lead 
to a spin period of about 1 ms after cooling and contraction. Assuming angular momentum 
conservation during contraction of the original proto-neutron star, that Ipns = kMR 2 PNS with 
k = const ~ |, and that Rns shrinks from 50 km to 10 km during the deleptonization phase, we 
take: 

'50 km\ 2 



p ™^lms^^^j -25 ms. (6) 

We note that this spin rate for the proto-neutron star implies a rather rapid rate of spin for 
the original iron core. Conservation of angular momentum (assuming conservation of mass and 
that the form factor for the moment of inertia is unchanged) gives 

Mp N s ( Rpns\ 2 nrM 

^Fe ^ -TT -5 Vp NS ~ 0.04s, (7) 

M Fe V K Fe J 

or a period of about 2.6 minutes. This represents a rotational energy of about 10 47 erg, an amount 
that is significant, but still much smaller than the binding energy of the iron core, ~ 10 51 erg. If 
the total helium core is in solid body rotation with the same period of about 3 minutes, then for a 
mass of 5 M Q and a radius of 2 x 10 10 cm, the rotational energy would be about 10 51 erg. This is 
clearly a demanding criterion. The specific parameters we adopt here in equation (g) and below 
are probably only reasonable if there is substantial differential rotation between the iron core and 
the surrounding star. 
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The moment of inertia of the proto-neutron star is Ipns — 3 x 10 46 c.g.s. With the fiducial 
period from equation (|6|) the rotational energy of the proto-neutron star is then: 

T-i 1 t 2 n m50 ( M NS \ ( ^PNS \ 2 ( RpNS\ 2 , c n 

Wns - 2 WW - 9 x 10 erg (^^j . (8) 

This is a substantial energy, but since only a fraction of it could be tapped it is not clear that this 
energy could power a supernova never mind a classical 7-ray burst. If the rotational period of the 
proto-neutron star is longer, this is even more true. The fraction of this rotational energy that can 
be tapped depends on the behavior of the magnetic field, and we turn to that subject. 



4.1. Magnetic Field Amplification 

The magnetic field of the proto-neutron star is uncertain. If the pre-collapse core has a field 
strength comparable to that of a magnetized white dwarf, ~ 10 6 G, then a field of ~ 10 12 G could 
arise from flux- freezing. For Ppns — 25 ms (equation ||) and r conv ~ 1 ms, the Rossby Number 
Rq ~ Ppns / T conv — 25 3> 1 is too large to allow an a — O type dynamo to operate. An alternative 
magnetic field amplification mechanism is linear amplification (Meier et al. 1976, Kluzniak & 
Ruderman 1998). In this process, the differentially rotating neutron star could wrap the poloidal 
seed field into strong toroidal fields which then emerge from the star through buoyancy. After 
revolutions of the neutron star, the initial seed (poloidal) field is wrapped and amplified to 
produce a toroidal field 

~ 2vrn S p , (9) 

where B p is the initial seed poloidal field. Buoyancy will operate to expel the field if the amplified 
final field Bf satisfies 

B 2 

- fBPcl (10) 

where fp — 0.01 is the fractional difference in density between the rising flux tube elements and 
the stellar material. Assuming sound speed c s ~ c/3 and density p ~ 10 13 g cm -3 , one finds 

B f ~ 2 x 10 16 G fl^pH 2 , (11) 

where Jb~2 = /b/0.01 and pi$ = p /10 13 g cm -3 . Assuming that the magnetic flux tube (a 
torus) occupies a volume of Vb/Vpns — 0-1 where Vpns — ~r ^Pivs 1 ' ^he energy contained in the 
magnetic flux tubes at the buoyancy limit is estimated to be: 

B 2 

E B ~ 0.1 x V PNS x -L (12) 

47T 

~ 1.6 x 10 51 erg. 

This magnetic energy is ejected from the neutron star in the rising magnetic flux tubes. 
This energy is comparable to the proto-neutron star rotation energy (equation ||). For the 
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adopted parameters, the proto-neutron star would have to lose substantial rotational energy to 
the magnetic field before the field would float on dynamical time scales. This magnetic energy is 
still likely to escape from the proto-neutron star, but the details may be complex and involve the 
subsequent contraction of the neutron star. 

The number of revolutions to reach Bf ~ 2 x 10 16 G is: 

"' = ff^ 3x103 (10^0)" < 13 > 

For Bq, ~ 10 12 G, the amplification time scale before the field is expelled by buoyancy is thus 

tf ~ n f Pp N s ~ 25 ms x 3 x 10 3 ~ 75 s, (14) 

at the beginning of the contraction of the proto-neutron star. As the neutron star contracts, it 
spins up. The timescale for the linear field amplification, tf, will decrease and the rotational 
energy of the neutron star will increase (cf §6). By the time the neutron star is spinning with a 
period of 1 ms, tf would be substantially less than the cooling, contraction time of ~ 10 s. Thus 
sometime during the contraction a dipole-type strong field is expected to be produced once the 
~ 10 16 G field emerges from the surface. A dipole field strength of ~ 10 14 G is expected from the 
random sum of flux tubes of ~ 10 16 G (Duncan & Thompson 1992; Thompson & Duncan 1993). 
It is likely that during the contraction phase a substantial fraction of E ro t t pNS — Eb,pns — 10 51 
erg would be extracted from the neutron star. 



4.2. Energy Transport 

We now consider the possible energy transport mechanisms from the neutron star to the outer 
stellar envelope. For P ~ 25 ms, the radius of the light cylinder is (cf. equation ||) Rlc,pns — 10 3 
km, comparable to the initial radius of the iron core and much larger than the proto-neutron 
star. For R < Rlcpns-, a rotating dipole field could exert a magnetic torque on the plasma 
around the star. Typical collapse calculations (Burrows, Hayes & Fryxell 1996; Mezzacappa, et al. 
1997) give for the density in the vicinity of the standing shock p ~ 10 8 g cm -3 with a pre-shock 
velocity of v ~ 10 8 cm s _1 . Using those as fiducial values, the characteristic Alfven radius for the 
proto-neutron star is (equation ^) : 

Ra * 150 km (ilm) {W^h) ' ( lOSg'cm-a ) ' ( l0» cm s~i ) ' ' (15) 



Equation ( |l~5l ) implies that for Bp^s substantially less than 10 14 G the Alfven radius would be 
less than the radius of the proto-neutron star and hence that the magnetic torque would be 
ineffective. For Bp^s ^ 10 14 G, the torque could have a significant effect. Equation ( |i~5"l) shows 
that Ra <C Rlc — 1000 km, so that a dipole field could be maintained over a significant volume. 



For the fiducial values chosen in equation (15), Ra is interestingly close to the stalled shock radius 
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~ 200 km, thus justifying the choice of density and velocity, to which Ra is not, in any case, very 
sensitive. 

There are two relevant timescales pertinent to the action of the torque. One is the time to 
deflect infalling matter from radial infall to substantially azimuthal flow. The second is the time 
for this torque to spin down the proto-neutron star and hence to deposit energy in the infalling 
matter. The time for the torque to deflect infalling matter is approximately the momentum per 
unit volume of the infalling matter times an appropriate lever arm, Ra, divided by the torque per 
until volume, ~ B Z (R)B^(R) ~ B 2 (R) (cf. Shapiro & Teukolsky 1983). This effect of the torque 
thus operates on a time scale 

pRav rRav 1 Ra n „x 
^ ~ 52 - Anpv 2 ~ 47r v > l ib J 



0.01s 



( RPNS \ ( B PNS \ V3 / p \ -1/6 / v N "4/3 

V5o k m y vio 14 g; vio 8 gcm- 3 y vio 8 



cm s~ 

This time is sufficiently short that the torque could substantially alter the flow of matter, 
preventing the radial infall interior to the shock that is common to all spherically symmetric 
models of core collapse. Multi-dimensional models that produce non-radial circulation flows in the 
matter beneath the standing shock would also be substantially affected for conditions similar to 
those reflected in equation (|i~6|). 

The time scale to extract energy from the proto-neutron star and deposit that energy in the 
surrounding plasma can be estimated by equating the torque, N, on the proto-neutron star to that 
on the gas, 

N = m~ [ R 2 B z {R)BJR)dR~ R A B{R A ) 2 . (17) 
JRa 

Defining the timescale appropriate to this action of the torque to be tt or ,2 = £l/Cl and using 
equation ( |i~7|) gives: 



1 SIpnsMnS / 10 x 

b^ir RpNsBpNsp 1 '^ 

170 s ( QpNs ) ( Mns ) ( RpNs \ ( ^ r 1 / n x ~' 2 ' 71 x ~' 



\250 s~ l J \1.5M Q J \50 km/ \10 14 GJ V10 8 gmcm-V Vl0 8 cm S - 1 
The rate of energy deposition in the plasma by this torque is thus: 

L tor = mn~2^nR PNS B PNS p 1/2 v, (19) 

fipws \ / Rpivs; \ 3 / Bp/vi \ ( \ 



1 x 10 49 ere s" 1 ( ^ PNS \ ( RpNS \ ( B P*s \ ( P 

ix 10 ergs { 25Qs -i){ 50km ) [ w u G ){ w s gi 



g cm" 3 / \10 8 cm s _1 
Note that since the torque is independent of the angular frequency will tend to decrease linearly. 
The timescale ttor,2 is sufficiently long and the associated energy deposition rate in equation ([l9|) 
is sufficiently small that it is unlikely that this energy deposition will directly contribute to any 
supernova explosion. On the other hand, the selective deposition of this energy very near the 
stalled shock might have a leveraging effect belied by the small deposition rate that could help to 
reinvigorate a stalled shock in conjunction with other effects such as neutrino deposition. 
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5. The Effect of an Axial Jet 

If, as in LeBlanc & Wilson (1970; see also Miiller & Hillebrandt 1979; Symbalisty 1984), 
an MHD jet is formed during the collapse phase, the maximum jet power is estimated to be 
E TO t/T(iyn — 10 51 erg/1 s ~ 10 51 erg s _1 . Although the details of the jet dynamics and energetics 
are not known, roughly ^ 10 51 erg s' 1 of power could be directed along the jet axis. The LeBlanc 
& Wilson calculation was criticized by Meier et al. (1976) as requiring extreme parameters of the 
progentor star. These issues need to be re-examined in the current context, but we note several 
things about the Meier et al. analysis. They argue that the MHD axial flow found by LeBlanc & 
Wilson will not propagate to the stellar surface as a jet. The calculation of Khokhlov et al. (1999) 
shows that this is not necessarily correct, at least for a helium core. Meier et al. based their 
analysis on stellar evolution calculations of the day, but they adopted a stellar core with central 
density of about 10 10 g cm' 3 giving a binding energy of about 10 52 ergs. This exaggerates the 
binding energy of the initial core by about a factor of 10 compared to modern calculations and 
gives an incorrectly small value of a key parameter of Meier et al., the ratio of the binding energy 
of the newly-formed neutron star to that of the initial core. Meier et al. also did not consider the 
possibility of an a — dynamo that could lead to exponential field growth (Duncan &: Thompson 
1992). Symbalisty found that both substantial rotation and magnetic field were necessary to 
affect the explosion and that the presence of the magnetic field and associated losses allowed 
deeper collapse. The whole question of the initiation of MHD jets in association with neutron star 
formation needs to be considered anew. Here we consider some basic properties of jet propagation. 

The jet would be stopped by the envelope when the jet is unable to provide the power to 
move the envelope material at a speed comparable to the jet velocity. This can be expressed as: 

L jet ~ R 2 An x p env v 2 jet x Vjet, (20) 



where Af2 is the solid angle of the jet. The jet would then be stopped in a length 

1/2 

Rst^ ' 



AQp env v 3 



(21) 



u jet_ 

In order to penetrate the star, the energy injected into the envelope at R & R s t should be enough 
to unbind the region of the outer stellar mantle impacted by the jet. The amount of stellar 
material impacted by the jet is 

AM env = M env — ~ 8 x W~ 3 M env f — J . (22) 

Unbinding this mass requires 

GM env M NS An 
L ^ tZ AnR env At ' (23) 



where At is the duration of injection of the jet. Equation (23) can be expressed as: 



3et S \IM Q J \1.5Mq J \W 5 kmj \0.1 J \ls J 



(24) 
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A jet of ~ 10 51 erg s 1 thus gives ample power to unbind a portion of the envelope occupying 
about 0.1 sterradian. 

The dynamics of the jet and impacted envelope material will depend on whether the mass of 
the jet is greater than or less than the mass of the displaced stellar envelope. The speed of the jet 
will be 

( 2L jet At \ 1/2 

V ^ 2c {-w-^) ■ (25) 

This implies a relativistic velocity for LAt ~ 10 51 ergs and Mj e t ^ 0.01M Q . If a comparable 
amount of energy is put into the displaced envelope material, the impacted envelope mass, AM env , 
would expand at a speed 

^"(^^) 1/2 - (26) 
If this matter were displaced radially out of the star, it could also acquire relativistic speeds for 
LAt ~ 10 51 ergs and AQ & 0.1. In practice, the displaced material will tend to be accelerated 
sideways by shocks induced by the passage of the jet and the energy will be shared by the whole 
envelope at roughly the sound speed in the envelope. 

The propagation of the jet requires more study, but the calculation of Khokhlov et al. (1999) 
gives some qualitative insight into the behavior to be expected. As the jet propagates, a bow 
shock runs ahead of it. The speed of the bow shock is less than that of the matter inflow into the 
jet, a ratio of about 0.5 to 0.7 for the particular case explored by Khokhlov et al., but that ratio 
will depend on the speed and density of the jet, its opening angle, and the structure of the star 
through which the jet propagates. 

The bow shock of the jet will both heat material and cause it to expand sideways. The 
opening half angle of the jet will then be approximately 

/i V env V env S Q V env S /r>i-r\ 

~ ~ 0.1 rad — ~ 5 — . (27) 

Vbow v bow,9 ^bow,9 

For Venv ~ 2 x 10 8 cm s" 1 and v bow ~ 2 x 10 9 cm s _1 , AO ~ 5° or AVL/Ait ~ 0.004. The actual 
dynamics of the jet will depend on nested cocoon-like shocks from the bow shock and subsequent 
expansion, as illustrated by Khokhlov et al. (1999). 



6. The Neutron Star Phase 

Concurrent with the propagation of any MHD jet and magnetic torque exerted by the 
proto-neutron star, the neutron star will contract and new physics can come into play. We assume 
that the cooling and contraction leads to a neutron star rotating at a period of about 1 millisecond. 
This implies a much smaller radius for the light cylinder, 

Rlc = ^~ * 50 km (p™) . (28) 
Q NS Vims/ 
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This is only a little bigger than the neutron star radius. If the field remains about the same as 
the proto-neutron star, B ~ 10 14 G, the Alfven radius will be comparable to Rlc-, but if the 
field amplifies during the contraction, as is plausible, the Alfven radius will expand to become 
substantially greater than Rlc- Under this circumstance, a stationary dipole configuration at and 
beyond the Alfven radius can no longer be maintained. 

During the contraction phase, the Rossby number decreases in proportion to the rotational 
period. Assuming the neutron star convective time scale remains at about 1 ms, the condition 
Ro 1 will be reached if the neutron star ends up with a period of about 1 ms. This means that an 
a — Q-type dynamo could be initiated that would build a strong magnetic field up to ~ 10 17 — 10 18 
G (Duncan & Thompson 1992; Thompson & Duncan 1993; equipartition with the rotation would 
give ~ 10 18 G). In this extreme circumstance, the surface field could be ~ 10 15 — 10 16 G for the 
dipole configuration near the neutron star surface (in the absence of the light cylinder). This 
phase of field growth occurs exponentially in contrast to the linear amplification associated with 
the relatively slowly-rotating initial proto-neutron star. If the dipole field grows to Bjys — 10 16 G 
as the radius shrinks to R^s — 10 6 cm, then, all else being equal, the Alfven radius will be ~ 140 
km, substantially larger than Rlc, as just noted. If the ram pressure has declined, Ra will be 
even larger. 

The rotational energy after the contraction is 

Era , s , lr NS n h , e x up erg (^-) (J^) 2 (J^) 2 . ( ») 

Setting aside the baryon-loading problem for the moment, this energy is comparable to the largest 
energy associated with any 7-ray burst, ~ 3 x 10 54 AO/4tt erg for GRB 990123 (Kulkarni, et al. 
1999) for a degree of collimation, AQ/4tt & 10 -2 . This degree of collimation has been deduced 
for some afterglows, in particular for GRB 990123 itself, and as noted in the previous section, is 
about the order expected for a matter-dominated jet (see also Khokhlov et al. 1999). The degree 
of collimation of the subsequent flow of electromagnetic radiation is unclear. We return to that 
topic below. 

The rotational energy of the contracted neutron star is radiated away in the form of intense 
waves of frequency £Ins generated at the speed of light circle. To generate traditional propagating 
electromagnetic waves as opposed to MHD waves, the displacement current must exceed the 
plasma current. This criterion can be written roughly as 



where e is the electron charge and iVo is Avogodro's number. The density within the star vastly 
exceeds this limit, so the energy will be generated as MHD waves. These waves will have extreme 
properties that require futher study. In particular, the Alfven speed can be estimated as 

VA = V ^ ~ 2 x 10 11 cm s- 1 ( J?™) ( P ,) . (31) 

\8717J vio 16 g; vio 8 g cm- 3 ; v ; 
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If the surface field is of order 10 G, then the field at the speed of light circle may be about 
10 15 G. Even for weaker fields, the density near the speed of light circle will decline in dynamic 
reaction to the large deposition of pulsar energy. Equation ([H]) thus serves to argue that this 
simple Newtonian expression for the speed of MHD waves is likely to be invalid since it is 
neither relativistic nor applicable in the strong field limit where the waves are not a first order 
perturbation to the plasma. Instead one must treat the generation and propagation of high 
amplitude, ultrarelativistic MHD waves (UMHDW), a task beyond the scope of this paper. 
Rather, we will sketch the basic energetics and the potential behavior of these waves. 

If the density gets sufficiently low as the UMHDW propagate or because of dynamical reaction 
to energy deposition that creates a cavity causing the density to drop, then the plasma current 
can be exceeded by the displacement current and the waves can be described as Large Amplitude 
Electromagnetic Waves (LAEMW; Usov 1994; Blackman & Yi 1998). The properties of these 
waves are also not completely understood. Any LAEMW cannot propagate through a plasma if 
the plasma frequency 

1/2 

~ 9 x 10 3 n^ 2 Hz, (32) 
exceeds the LAEMW frequency Qns- The corresponding condition on the electron density is: 

^ lcm-3 (loT^r) 2 - ^ 

For densities exceeding this value, the LAEMW would have a very small skin depth and cannot 
propagate in the plasma. They would effectively be reflected by the plasma. The density exceeds 
this critical value under any stellar conditions. 

In the discussion below we neglect the production of traditional Alfven waves as an energy 
loss mechanism. Alfven waves occur when the magnetic field is relatively weak and the waves 
are a propagating perturbation. Alfven waves will thus be generated beyond the magnetosphere 
where the magnetic field no longer dominates the particle dynamics. In the proto-neutron star 
case, Alfven waves could be generated in the region beyond the magnetosphere and within the 
extended light cylinder. The torque acting at Ra probably dominates those losses, but this is 
worth more detailed consideration. Once the neutron star contracts and the light cylinder passes 
within the magnetosphere, where Alfven waves are not a relevant concept, the dominant losses 
will occur at the speed of light circle by generation of UMHDW. Alfven waves might be generated 
up the rotation axis within the light cylinder but beyond the Alfven radius. The global behavior 
of the open magnetic field lines that cross the light cylinder in the complex environment we are 
discussing is unclear. The dynamics of these open field lines again might lead to the production of 
Alfven waves and associated losses, but it seems likely that in all these possibilities, Alfven waves 
remain a secondary process. 
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The luminosity in UMHDW is estimated to be 

47T K LC C 

assuming the UMHDW to be generated at Rlc and the magnetic moment of the neutron star 
to be /xjvs' = BnsR%s- -F° r t ne conditions of the contracted neutron star which has initiated an 
a — f2 dynamo, we expect 

T A 1052 -1 (RNS \ 6 (ENS VfJlNS V 

W^4xl0 ergs ( J ^j (^-^j , (35) 

which will last for a duration of 

E r ot,NS ^ 2 s ( — f ^ jvs ~\ ^ f ^NS \ 2 f ^NS 



tMHD " W " 2 S ll^o ) W^J I^gJ llO^J • (36) 

The UMHDW are likely to initially be bottled up near the site of their production near Rlc 
since the surrounding plasma can respond only at the subrelativistic sound speed. The UHMDW 
would push the plasma aside and seek the easiest way out. An obvious possibility is that they will 
"burn" a channel along the rotation axis, following the path of the previous MHD jet. We assume 
that the most intense production of UMHDW occurs after the completion of the contraction 
of the neutron star, and hence that it is delayed by the cooling time ~ 5 to 10 s, about the 
time necessary for the initial matter jet to propagate through the stellar core. A mass of about 
AM ~ M env ^- ~ 1CP 2 M env will be pushed aside by the initial matter jet. Whether the region 
impacted by the matter jet will be rarefied depends on the mass of the jet. If the jet is rather 
massive, as in the case of the calculation of Khokhlov et al. (1999), then the jet remains denser 
than the stellar environment as it propagates into the envelope. In this case, the jet acts as a 
"plug" until it disperses after several dynamical times. The UHMDW generated by the pulsar will 
tend to propagate through the lowest density regions. If there is a density minimum along the 
rotation axis, there will be a tendency for energy to flow in that direction. The dynamics are likely 
to be complicated and we will just sketch the possibilities. 

The UMHDW should be rapidly isotropized in the co-moving frame, so they can be considered 
as a relativistic gas while they are trapped within the stellar core. It is not so clear that they are 
thermalized since their characteristic wavelength 

Xmhd ^ cP NS ^ 2ttR lc ~ 300 km, (37) 

is comparable to or larger than the scale of the region in which they are generated. This is an 
important issue since if the energy is thermalized then the effective temperature will be, 

Tmho ^ (|^f ) V4 - 2 x 10" K E^Rf'\ (38) 
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where R3 is the radius in units of 10 3 km, so that much of the energy would go into the formation 
of pairs which might be dissipated by adiabatic expansion and difficult to recover. We will ignore 
this possibility for the moment, assuming that the energy cannot thermalize. 

If the UHMDW act like a relativistic gas, then they will exert a pressure of 

Pmhd - \ (|^f ) - 8 x 10 26 erg cm" 3 E 52 R^. (39) 

This may be contrasted to the pressure required to ensure hydrostatic equilibrium in the gravity 
of the neutron star which we can represent crudely as, 

Phse - - 1 x 10 28 erg cm" 3 ( ^f] ' R^. (40) 



R 4 & \1.5M . 

By this measure, the pressure of the radiation will be less than required to produce HSE for small 
radii, and will exceed the pressure corresponding to HSE at a radius of 

Rhse * 1 x 10 4 km {J^jEg. (41) 

The regime where the UMHDW begin to dominate will, of course, depend on the density profile in 
the collapsing matter and the whole environment is, in any case, dynamic so that considerations 
of HSE give only a partial perspective. 

The production of the UMHDW near the light cylinder may lead to strong Rayleigh- Taylor 
instability as the relativistic fluid pushes on the surrounding gas. The subsequent behavior will 
depend on the rate at which the UMHDW "gas" expands quasi-spherically in Rayleigh- Taylor 
fingers compared to the rate at which the waves will selectively propagate upward along the 
rotation axis. It seems plausible that a substantial majority of the energy will flow up the axis, 
the path of least resistance. Another possibility is that the energy density becomes so high 
compared to the surrounding matter, near the neutron star or further out in the mantle (equation 
41), that the stellar matter becomes an insubstantial barrier to the nearly free expansion of the 



UMHDW bubble. Since the total energy in the UMHDW could be as high as 10 52 erg, this is 
a real possibility. Whether the UMHDW remain collimated may thus depend sensitively on the 
timescale on which they are generated and hence the instantaneous energy density throughout the 
star. These possibilities clearly need to be investigated with appropriate numerical calculations. 

The result of the generation of the UMHDW could be rather different if they are thermalized. 
In this case, the UMHDW would be replaced with copious pairs (cf. equation ^). The pair 
gas would also act like a relativistic gas with the same issues of dynamics just discussed. In 
addition, there would be the complications of annihilation at boundaries with normal matter and 
the possibility of strong adiabatic losses after a pair jet broke out of the star and underwent free 
expansion. Such loss of pair thermal energy to rest mass and kinetic energy might be at least 
partially recovered if the pairs interacted with a surrounding wind of normal matter (see §8). 
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7. Evolution of Energy Emission and Gamma- Ray Bursts 



There are two stages when 7-rays might be generated. The first is when the bow shock that 
proceeds the initial mass-dominated jet impacts on the stellar photosphere. The second phase is 
when a collimated flow of UMHDW erupts from the surface at about the same time. We consider 
those in turn. 



Some 7-ray emission may be generated by the phase of shock breakout as the shock associated 
with the initial matter-dominated jet runs down the stellar density gradient and breaks through 
the photosphere. An associated mechanism is for the accelerated matter to reach relativistic 
speeds and then collide with some external matter. 

As the bow shock that proceeds the jet runs down the exponential stellar density gradient in 
the photosphere it will accelerate, in turn heating and accelerating the matter, a mechanism to 
produce 7-rays first described by Colgate (1974, 1975; see Matzner &; McKee 1999 for an updated 
discussion of analytic models of shock breakout). The question of how this matter expands and 
radiates and perhaps collides with an external environment is beyond the scope of this paper, but 
worth more detailed study. Here we remark on the basic properties at the time of shock passage. 

Energy is deposited in the atmosphere by the bow shock. The bow shock, in turn, derives its 
energy from the jet. The power carried by the jet is ~ l/2pj e tVj et AA, where AA is the fractional 
area of the star that is impacted by the bow shock. This power will be delivered to the photosphere 
in a time ~ l p hot/ v jet, where l p hot is the depth of the photosphere. The energy deposited in the 
photosphere during the bow shock break out phase is thus: 



where Vjet,io is the velocity of the bow shock in units of 10 cm s , l p hot,7 is in units of 10 cm 
and is in units of 10 18 cm 2 . For the calculation of Khokhlov et al. (1999) the jet density is 

about pj et ~ 10 3 g cm -3 , and the velocity steepens to about 90,000 km s _1 as the jet approaches 
the outer density gradient of the helium core at which point the resolution of the density profile 
degrades. For pj et ~ 10 3 g cm -3 , Vj e t t io ~ 1, Iphotj — 1, and AAi$ ~ 1, the energy would be 
~ 5 x 10 ergs, comparable to the 7-ray energy in the burst of SN 1998bw/GRB 980425. 

The corresponding temperature, making the crude, and not necessarily correct, assumption 
of thermalization to a radiation-dominated gas, is 



7.1. Bow Shock Gamma- Ray Emission 




j et l phot AA ~ 5 x 10 45 erg p jet v 2 jet l0 l phot jAA 1& 



(42) 




(43) 



For p jet 
7-rays. 



10 3 g cm 3 , Vjet.io — 1) equation fl43"D gives T ~ 10 10 K with the possible emission of 
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This hot, accelerated matter would then expand and radiate. Whether it could account for 
the 7-ray burst in SN 1998bw/GRB 980425 will require more careful consideration. We note that 
if this energy accumulates in the density gradient of the photosphere the matter will be optically 
thin, so not susceptible to adiabatic losses while it undergoes free expansion. The photosphere is 
already optically thin by definition and the opacity will be decreased by Klein-Nishina corrections. 
On the other hand, the jet is denser and hence more opaque than the photosphere so it is important 
to know where the energy represented by equation (|42|) resides. Although the conditions could 
vary widely depending on the nature and propagation of the jet, for the calculation of Khokhlov 
et al. (1999), the bow shock material is matter ablated from the jet and the characteristic density 
is about the same as the jet, p^ow — 10 3 g cm -3 . 

The emission properties of the jet/bow shock/stellar atmosphere region as the jet impacts the 
atmosphere thus require more careful study. Nevertheless, considerable hard radiation could be 
emitted before the phase of homologous expansion is reached. The time scale for this energy to be 
radiated is also of importance. The shock break out time, ~ 10 -3 l p hotj/ v jet,io s, is far too short 
to correspond to an observed 7-ray burst, specifically that in GRB 980425, but the controlling 
time scale will be the radiative timescale. 

Another possible means of producing 7-rays from the matter jet is to accelerate the matter in 
the photosphere to relativistic speeds and for it then to subsequently collide with a surrounding 
medium, perhaps a stellar wind. The mass fraction that can be accelerated in this way for an 
explosion with an impulsive energy input that drives a single shock through the outer stellar layers 
has been evaluated for CO and He stellar cores by Woosely, Eastman & Schmidt (1998). For their 
models, they establish a relationship 



where the parameter Q is independent of Lagrangian mass and radius in the ejecta, M e:r 32 is the 

mass external to the layer with Lorentz factor T in units of 10 32 gm, and we have taken (3 ~ 1. 

Woosely et al. find that Q ~ 3 x 10 5 . At the risk of overinterpreting their results, we note that for 

two otherwise identical models that differ only in the explosion energy input (models C06A and 

1/2 

C06C), Q scales with the explosion energy, E ex , approximately as E ex . In the following we have 
adopted Q ~ 10 5 £^ 2 51 , where E ex ^\ is in units of 10 51 erg. 

With this scaling, we find the kinetic energy in the homologously expanding material (e.g. 
after the immediate shock heating phase discussed above) to be 



This expression is for a spherical explosion. It is clear that the energy in matter accelerated to 
relativistic speeds is insufficient to account even for the weak 7-ray burst imputed to GRB 990425, 
as concluded by Woosely, Eastman & Schmidt (1998). 




(44) 



KE = TM ex c 2 ~ 1.6 x 10 44 erg T 



-5.66 P2.83 



(45) 



In the present context, it is relevant to estimate the difference if the shock did not propagate 
spherically, but were collimated. In this case, we are interested in the kinetic energy in a jet with 
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solid angle AQ/4tt or KEj\n = KEAQ/4tt. The relevant input energy is the energy in the jet, 
Ejet = (AQ/4:Tr)E ex , where E ex is now to be thought of as the equivalent isotropic energy of the 
collimated jet. Making this substitution in equation ((45|), we find: 

KE An ~ 7.2 x 10 47 erg T' 5 ™ f-jj^Ef^, (46) 

where f C oii-2 = AQ,/4ir in units of 10~ 2 . Equation fl46| ) shows that even with substantial 
collimation, the amount of energy put into mildly relativistic matter (r <; a few) by this 
mechanism is still insufficient to account for the putative 7-ray burst in SN 1998bw for an energy 
typical of the matter jet we have discussed here, Ej e t 51 — 1. There might be sufficient energy for 
Ejet,5i — 10, but even this energy would be insufficient to boost enough stellar matter to T ~ 10 
to account for SN 1998bw. 

From this analysis, we conclude that, even if it is strongly collimated, the physical process of 
running a single impulsively-induced shock (e.g. Colgate 1975; Matzner & McKee 1999) down the 
photospheric density gradient of a compact stellar core, accelerating that matter to relativistic 
speeds, and slamming it into a surrounding medium is unlikely to generate a substantial 7-ray 
burst of the kind associated with SN 1998bw. To produce a cosmic 7-ray burst with an energy of 
10 52 erg (isotropic equivalent of 10 54 erg) and T ^ 100 is out of the question. If jets from stellar 
collapse to produce either neutron stars or black holes are to generate 7-ray bursts in SN 1998bw, 
never mind the high redshift events, then the physics of the 7-ray production must be associated 
with the initial shock breakout phase (as discussed above) or by the prolonged emission phase of 
the jet passing through and beyond the photosphere. 



7.2. The Contraction/UMHDW Phase 

We now consider the time evolution of the contracting neutron star and the associated energy 
emission mechanisms implied by this scenario. 

When Rlc ^> RAi the magnetic torque is the dominant energy/momentum transport 
mechanism (§4.2). When Rlc *C Ra, Poynting flux/UMHDW transports the energy/momentum. 
In our scenario, the proto-neutron star continues spinning up due to contraction even though it 
transfers angular momentum through magnetic torque. Initially, when Rlc ^ Ra, the rate of loss 
of energy to electromagnetic Poynting flux would be (cf. equation |35|) 



t n m46 -if RPNS\ 6 ( BpNS \ 2 ( ^PNS \ 4 

— 2 x 10 erg s j i^-^j \2l30~lFv ' (47) 

This will be small compared to the energy deposited by the torque of the spinning neutron star 



acting at the Alfven radius (equation IS), 



! ln 49 -1 f VPNS \ f RpNs Y f B PNS \ ( p \ 1/2 ( V \ 

-— 1Xl ° 6rgS 1 250^ J UknJ U^A lOSgcm-a j [ W cm s^ J ' (4§) 
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After contraction when the condition Rlc ^ R-A is reached, we have: 



Rns V10 16 G/ V10 8 gcm- 3 y V10 8 cms- 1 

and 



and hence 



^£ ~ n 9 f fijvs ^ ~' f ^ s ^ _1/3 f p - \ V6 f ^ ^ 1/3 ( Uns 1 i - 1 1 



where p and u are to be evaluated at the Alfven radius. Since energy injection during the phase 
when Rlc > Ra will drive the density down in the vicinity of the magnetopause and the standing 
shock, and Bn$ is likely to increase from 10 14 G to 10 16 G continuously (but suddenly) during 
the contraction by the linear amplification mechanism and dynamo, the epoch during contraction 
when Rlc — Ra will occur is difficult to estimate; however, since essentially all the energy 
emitted will be used to expel the envelope (i.e. either by torque or UMHDW) a more significant 
explosion/expansion is expected when Rlc & Ra- At this epoch we expect that if the torque were 
still active the rate of loss of rotational energy due to this mechansim would be 

,49 .„ .-i ( Vns \ ( Rns \ 3 / B NS \ / p \ 1 ' 2 



L to ^3xl0 ergs- ( J ^j j ( 10 8 g cm -3 j [ w s cm s -i j • ( 52 ) 

This power would in any case be overwhelmed by the loss of energy to the UMHDW with the 
enhanced rotation and magnetic field, 



This power output will last for a few seconds (equation 36). The injected energy is then <; 10 52 
erg, which is enough to drive a significant shock. In this scenario, both a 7-ray burst and an 
asymmetric supernova could be produced during the "magnetar" phase. 

The UMHDW will likely flow up the rotation axis starting with a cross-sectional area of 
~ R\q and rising as an intense photon "bubble." Assuming that this bubble pushes matter aside 
at roughly the local sound speed, the opening angle will be: 

Omhd ~ — ~ 3 x 1(T 3 rad ~ 0.2°, (54) 
c 

where c s is the sound speed in the helium core, ~ 10 8 cm s . This corresponds to a solid angle of 
AJ1/47T ~ 2 x 1CU 6 . The UMHDW may then propagate up the axis in a radiation-dominated jet 
with smaller cross section than the original matter-dominated jet. If this is the case, the channel 
carved by the UMHDW may remain substantially smaller than the characteristic wavelength of 
the UMHDW until the waves break out of the surface of the star. 
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To avoid the baryon loading problem and to generate a 7-ray burst with large Lorentz factor, 
the hole through which the bulk of the UMHDW emerge should contain less than 

AM & (55) 

(Rees & Meszaros 1992). For T ^ 100, the requirement is AM <, 10 27 g ~ 5 x 10~ 7 M . It is 
difficult to estimate whether the UMHDW jet will entrain such a small amount of matter. There 
are two qualitative possibilites for the propagation of the UMHDW jet. There may be a density 
minimum along the axis of the previous matter jet as suggested by the calculations of Khokhlov et 
al (1999). If the UMHDW propagate through such a low density axial channel, they will naturally 
emerge somewhat more collimated than the original jet. On the other hand, if the jet is denser 
than the surrounding stellar matter and the jet acts like a plug on the axis, the UMHDW might 
propagate in a cylindrical blanket around the plug. They would then emerge with an annular cross 
section. This would result in yet another complication for predicting the observational aspects of 
such an event. 

If the matter jet leaves behind a region of very low baryon density, or the UHMDW emerge 
beyond the end of the more slowly propagating matter jet, then the propagating UMHDW could 
enter a low density region in which the density of the environment is lower then the critical 
Goldreich-Julian density (Goldreich &; Julian 1969; Shapiro & Teukolsky 1983) of 

PGJ * 10- g cm- ( -») (t^V 1 , (56) 



V 10 16 G J V 1 ms 



where the flux- freezing and force free conditions are broken (see equation 30). At this point, 
the transition from UMHDW to LAEMW can be made. Reconnection of magnetic fields and 
rapid acceleration of particles in a pair plasma follows (e.g. Asseo, Kennel & Pellat 1978; Usov 
1992; Michel 1984; Thompson 1994; Blackman, Yi & Field 1995). The pair plasma could in 
principle be accelerated to a high bulk Lorentz factor up to ~ 10 5 , although the exact value of the 
maximum bulk Lorentz factor could be significantly lower than this due to baryon loading and 
other complicated factors reducing acceleration efficiency (cf. Usov 1992, 1994). 

The rest-frame 7-ray luminosity would be determined by the electromagnetic power, i.e. 

* aw * 4 x ^ erg (A) ( y (_^)« ( jw) 2 , (57) 

where / 7 is the 7-ray emission efficiency factor. During this phase, the pulsar rotation energy is 
used to power the 7-ray burst while the spin of the pulsar slows down as 

,4 -if 



VnS = IO^-Mt^t) (58) 
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where Qns,i is the initial spin frequency of the neutron star and i e m,i is the initial time of the 
spin-down phase driven by the electromagnetic dipole radiation emission. The time evolution of 
the 7-ray luminosity is then given by 

t a ir>51 -1 ( fl \ ( ^NS,i \ 4 / RNS \ 6 / B NS \ 2 /cn x 

L ° bs * 4X10 6rgS [u)[W^J UokmJ [W^gJ (59) 

1 + 6 x icr 2 ^ Mns ^ Ns ' i ^ 2 ( ^ ns ^ 4 ( Bns ^ 2 f t ~ te 



i.5M ; vio 4 s- i y viokm; vio 16 G7 V is 

The LAEMW-powered 7-ray emission stage shows a characteristic initial phase during which the 
luminosity remains nearly constant. This phase is expected to be followed by a phase in which the 
luminosity decreases as oc t~ 2 . As shown in Blackman and Yi (1998), the 7-ray emission via the 
synchrotron-Compton process gives the luminosity of the peak oc We note that if the proton 
fraction changes after some complicated acceleration processes, then the synchrotron-Compton 
emission process and the above scalings could be significantly different. 

To get 7-rays directly from the LAEMW, there must be a density in the environment that 
is below the Goldreich-Julian density so the currents cannot be supported in the plasma and 
a pair cloud is spontaneously generated. The Goldreich-Julian density scales with the ambient 
magnetic field. If that field falls off like in the propagating Poynting flux, then beyond the 
helium core at 10 6 km, the Goldreich-Julian density will be (equation |56| ) p ^ 10~ n g cm~ 3 . 
Densities this low might occur immediately beyond the star (and the jet) if the star is embedded 
only in the interstellar medium, but it is likely that the progenitor of a Type Ib/c supernova is 
surrounded by a wind or other mass resulting from mass-loss processes. For a constant velocity 
wind at 10 8 t?g cm s -1 carrying mass at a rate 10 -5 M_5 M yr _1 , the density is 

p = 5 x 1(T 9 g cm" 3 M^v^R^r 2 , (60) 

with radius in units of 10 5 km. This density will be less than the Goldreich-Julian density for 
10 6 km, several times the radius of the helium core. 

Although the possibility of 7-ray emission by pair cascade is not ruled out, it is important 
to consider how the UMHDW/LAEMW escape from the stellar core and associated processes. 
A sub-relativistic matter-dominated MHD jet leaves behind a baryon-rich environment. The 
common statement is that if the baryon loading is too high a 7-ray burst of observed properties 
cannot be produced. We show below that this is not necessarily the case. 

Since the plasma density (baryons, leptons, and photons) is high, the perfect MHD condition 
is always maintained throughout the phase when the UMHDW propagate within the stellar core. 
Due to the high radiation density and the large amount of high temperature stellar material, the 
optical depth for outgoing radiation is exceedingly high. This makes the 7-ray emission efficiency 
very low. Under these circumstances, the bulk of the energy supplied by the UMHDW may be 
used to accelerate shocks parallel to the rotation axis and the axis of the matter jet. After the 
shock breaks out, 7-ray emission could occur as the kinetic energy of the shock is dissipated into 
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particle acceleration and subsequent synchrotron-Compton emission. We note one important 
difference between this situation and that of the original Colgate mechanism. As noted in §7.1, 
in the Colgate mechanism there is an impulsive input of energy in the stellar core and a single 
shock that propagates outward. This leads to a single, short pulse of hard emission. This process 
may, indeed, work when the shock driven by the UMHDW first encounters the steep density 
gradient at the stellar surface or at the tip of the jet, but in the current situation the pulsar 
and the UMHDW/LAEMW continue to drive shocks for an extended time. The density profile 
will be altered as the continuing shock energy is deposited and that reaction must be considered 
self-consistently. 

It is difficult to estimate the 7-ray emission efficiency in this process. Even for a very low 
efficiency 55 10 -4 , however, a 7-ray burst event such as SN 1998bw/GRB980425 could be generated 
provided that the UMHDW power is large enough, <^ 10 51 erg s _1 as has been assumed. If the 
efficiency of production of 7-rays approaches unity, then, with appropriate collimation, a 7-ray 
burst detectable at cosmological distances could be produced. Aside from the complex details of 
the 7-ray production, the overall bolometric luminosity evolution is expected to follow the time 
evolution described by equation ( 5£ ) . 



In this context, we take note of circumstances when the deposition of energy into baryons 
does not necessarily doom the production of a 7-ray burst. As noted by Protheroe & Bednarek 
(1999), if protons can be accelerated to energies where pion production is efficient, then 7-rays 
can be produced in the subsequent pion decay. To reach this regime, the protons must be given 
at least mildly relativistic energies, in excess of 1 GeV. This requires that the energy of the 
UMHDW/LAEMW be shared with less than 0.005 M Q E 52 of baryons. The shock driven by the 
UMHDW/LAEMW should boost some protons to even higher bulk velocities and the energy 
deposited by the UMHDW/LAEMW will plausibly go into a mass less than the total mass of the 
precursor subrelativistic jet if their effects are concentrated in the relatively low density matter 
along the axis of the jet or the matter immediately surrounding the jet. Another alternative is that 
the protons will be accelerated by the momentum associated with the Poynting flux. Deposition 
of the energy of the pulsar spindown via UMHDW selectively into the bulk motion of protons 
with energies substantially above the pion production threshhold is thus a distinct possibility. 
Thus rather than being a handicap, the "baryon loading" of the jet could be an advantage in the 
conversion of pulsar energy into 7-ray energy. 

Although a quantitative calculation is required, we envisage a process by which the protons 
are accelerated to high energy with the associated efficient production of pions. To produce the 
pions, the protons must collide with a "target." In the present context, a good candidate for 
the target is the wind expected to be present around the progenitor. The stopping length for 
proton-proton interaction to produce pions is about 100 cm 2 gm" 1 . For the density given by 
equation ( |60[ ) , the column depth of the wind is 

1 = 50 g cm~ 2 M^ 5 Vs l Rs l , (61) 
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with radius again in units of 10 5 km. Thus a substantial wind might provide the stopping 
medium for the conversion of high-velocity protons to pions in the vicinity of the surface of a 
hydrogen-stripped star. 

The pions produced when the high-energy protons collide with the wind then decay and 
produce very energetic 7-rays. Note that the stopping length for 7-rays is about 30 cm 2 gm _1 , so 
a "target" that stops protons will be somewhat optically thick to 7-rays. These pion-decay high 
energy 7-rays do not immediately correspond to the observed 7-rays in a 7-ray burst. Rather, 
they induce a pair cascade through photon-photon collision. This pair fireball could then produce 
the observed 7-ray burst. The strength of a 7-ray burst produced in this way will depend on the 
relative efficiency of production of pions versus neutrinos. The latter would be an energy sink for 
this process. 

Note that equation ( |6l[) neglects any affect of the ambient or stellar magnetic fields or those 
associated with the UMHDW/LAEMW. If the field were strong and the interparticle spacing 
small, the shock generated by the UMHDW/LAEMW could be collisionless and there might be 
no pion production. Alternatively, if the shock accelerates protons into a small magnetic field 
then this estimate would be valid. As an example, if the field scales as R~ 2 from a central dipole 
of 10 16 G, then the Larmor radius near the surface of the star would be about 10 -4 cm. If the 
field is associated with the LAEMW and scales as then the Larmor radius would be about 
10~ 8 cm. At the density given in the wind by equation (|60|), the interparticle spacing is about 
10~ 5 cm. Thus gyrating protons might still collide in the 10 8 G dipole field with an effective 
path length enhanced by the gyration, and equation (^) could be a lower limit, or they could be 
rendered collisionless in the 10 12 G field associated with the LAEMW. On the other hand, since 
the UMHDW/LAEMW could dominate the dynamics of the protons, the idea of a shock among 
particles carried along by the LAEMW might be moot anyway (Usov 1999). Clearly, the physics 
associated with these processes is complex and in need of deeper understanding. The only point 
to be made here is that the physics of proton acceleration and pion production in the wind might 
be relevant and needs further study. 



8. Disscussion and Conclusions 

We have shown that the contraction phase of a proto-neutron star could result in a substantial 
change in the physical properties of the environment. The following parameters are relevant: 

• Ppns — 25ms • Pns — ^ ms 

• Rpns — 50 km • Rns —10 km 

• Ra,pns ~ 150 km • Ra,ns ~150 km 

• Rlc — 10 3 km • R LC —50 km 

• I PNS = kR 2 PNS M NS ~ 3 x 10 46 c.g.s. • I NS = kR 2 NS M NS ~ 10 45 c.g.s. 
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• R ~ PPNS/Tconv > 1 



— > ol — £7 dynamo 
• exponential growth of field 



— ► no dynamo 
• magnetic field amplified within 



~10s by linear amplification or 



flux-freezing 

• MHD jet & torque action 

• hole punched & envelope torqued 



• UMHDW generated: S = ^{E x B) 



• relativistic jet 

envelope accelerated and heated by waves 



When the rotating magnetized neutron star first forms there is likely to be linear amplification 
of the magnetic field and the creation of a matter-dominated jet, perhaps catalyzed by 
subrelativistic MHD effects, up the rotation axis. The energy of the proto-neutron star is sufficient 
to power a significant matter jet, but unlikely to generate a strong 7-ray burst. The matter jet 
could generate a smaller 7-ray burst as seems to be associated with SN 1998bw and GRB 980425 
by the Colgate mechanism as it emerges and drives a shock down the stellar density gradient 
in the absence of a hydrogen envelope, e.g., in a Type Ib/c supernova. As the neutron star 
cools, contracts, and speeds up, two significant things happen. One is that the rotational energy 
increases. The energy becomes significantly larger than required to produce a supernova and 
sufficient, in principle, to drive a cosmic 7-ray burst if the collimation is tight enough and losses 
are small enough. In addition, the light cylinder contracts significantly, so that a stationary dipole 
field cannot form and the emission of strong UMHDW occur. Tight collimation of the original 
matter jet and of the subsequent flow of UMHDW in a radiation-dominated jet is expected. 

The UMHDW will propagate as intense low frequency, long wavelength Poynting flux. They 
may be isotropized to act as a relativistic fluid, but not thermalized since they have a frequency 
much less than the plasma frequency. The UMHDW "bubble" could be strongly Rayleigh-Tayor 
unstable, but still may propagate selectively with small opening angle up the rotation axis as 
an UMHDW jet. Alternatively, the impulsive production of UMHDW could render the stellar 
matter nearly irrelevant as a confining medium. If a UMHDW jet forms, it can drive shocks which 
may selectively propagate down the axis of the initial matter jet or around the perimeter of the 
matter jet. The shocks associated with the UMHDW jet could generate 7-rays by the Colgate 
mechanism as they propagate down the density gradient at the tip of the jet or there could be 
bulk acceleration of protons to above the pion production threshold. The protons could produce 
copious pions upon collision with the surrounding wind, thus triggering a cascade of high energy 
7-rays, pairs, and lower-energy 7-rays in an observable 7-ray burst. Yet another alternative is that 
the UMHDW could eventually propagate into such a low density enviroment that they directly 
induce LAEMW and pair cascade (see also Thompson & Madau 1999). 
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The radiation-dominated jet cannot form for several seconds as the neutron star contracts, 
spins up, and generates a large magnetic field, but then it propagates faster than the matter- 
dominated jet. In this circumstance, the matter jet could precede or follow the radiation-dominated 
jet. In the former case an X-ray precursor could be generated. In the latter case the matter jet 
might not be conspicuous at all. 

There are a number of reasons why the processes we have outlined may not be as effective as 
we have assumed. One is that the rotation of the neutron star may prevent contraction to the 
high densities and rotation rates on the timescales we have assumed (Fryer & Heger 1999). This 
could affect the convection in the neutron star and hence the generation of the magnetic field 
by the a — 0, dynamo mechanism. This is a complex issue, of course, since the presence of the 
magnetic field will lead to energy loss and conditions of greater contraction, as we have invoked 
here (Symbalisty 1984). 

Another complication we have ignored is that the mean dipole field that forms may have its 
axis tilted with respect to the spin axis. This may not be a critical factor, since the subsequent 
dynamics may be dominated by the density distribution surrounding the neutron star that is 
predominantly set by the angular momentum, not the magnetic field. The question of what 
fraction of the pulsar energy goes to drive quasi-spherical expansion and what fraction propagates 
as co-linear UMHDW clearly requires greater study. We also noted that if the UMHDW are 
thermalized in the stellar core the result will be the copious production of electron/positron pairs. 
At first, such a pair cloud will behave like a relativistic fluid so there may be little difference in 
the dynamics. As Rayleigh- Taylor instabilities ensue, the pairs might get mixed with ordinary 
matter and the positrons annihilate. Even if the pair bubble escapes up the axis, as we envisage 
for non-thermalized UMHDW, it will expand as it propagates, especially after breaking through 
the stellar surface. Much of the thermal energy could then be lost to kinetic energy by adiabatic 
expansion. If there is no "working surface" with which this pair cloud could collide, then much 
of the energy could be lost. On the other hand, the sort of wind expected (equation |(]) would 
give a stopping length (equation ^l|) that will easily stop the pairs. The issue would then be the 
efficiency of conversion of their kinetic and rest mass energy into 7-rays. 

There is an interesting question of the opposite sort concerning the possibility that the 
production of large energy in UMHDW could be too effective. For instance, if an UMHDW jet 
could propagate through the naked core of a Type Ic, but would get stopped and dissipated in a 
red supergiant, then we might find that the explosion energy of Type II is systematically higher 
than that of Type Ic. The fact that there is no clear evidence for this may suggest that the 
production of 10 52 ergs in UMHDW is not a common occurrence in core collapse, but this question 
requires further study. 

The question of whether or not pulsar spin-down will produce a 7-ray burst depends on such 
factors as the initial rotation rate, the strength of the dipole field that evolves, the tilt of such 
a field compared to the spin axis, and the density of the progenitor wind. Issues of uncertain 
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physics aside, it is clear that this mechanism might not be robust in the production of 7-ray 
bursts, but might produce 7-ray bursts of varying strength depending on natural variation in the 
circumstances of a given collapse event. 

Any 7-rays emitted by any of these processs are likely to be strongly collimated. The 
luminosity of the emitted radiation will depend on the geometry of that emission. We have noted 
here that the energy produced by the spin-down of the pulsar could emerge from the stellar surface 
along the axis of a low-density matter jet, or in an annulus surrounding a high density jet. Either 
of these cases will give a Lorentz factor that depends strongly on the aspect angle of the observer. 
Computation of the resulting luminosity is thus distinctly non-trivial. 

Simple 7-ray burst models invoking collimation assume that the energy in a collimated burst 
scales simply as AH/Att compared to that deduced in an isotropic geometry. This assumes that 
the collimated jet nevertheless expands as the section of a sphere with a cross section scaling as 
R 2 . We note that if the relativistic flow is truly collimated, this may yet be an overestimate of the 
energy required to power a given observed 7-ray burst. The calculation of a jet emerging from a 
stellar core by Khokhlov et al. (1999) shows that the dynamical jet that precedes any possible 
7-ray burst is nearly linearly collimated and does not expand in cross section as R 2 . In the absence 
of an understanding of the dynamics of the subsequent flow of UMHDW/LAEMW, we do not 
know the geometry of the relativistic flow, but if the cross sectional area grows less strongly than 
R 2 , then the solid angle AQ/4tt will be a decreasing function of distance. The energetics (and 
hence the derived rates of occurence) of 7-ray bursts will depend not only on the fact, but the 
geometry of the collimation. If the cross section expands less rapidly than R 2 , then substantially 
less energy may be required to produce a given observed 7-ray burst. 

The observed nature of any 7-ray burst will depend on whether or not one is directly 
witnessing a strongly relativistic bulk flow. There is general agreement that the observed 
afterglows of the cosmic 7-ray bursts represent relativistic blast waves. The question of whether 
that is true or not for the primary 7-ray burst is still in contention with models based on internal 
shocks in a "central machine," in which the burst duration is the fundamental physical time scale 
of the energy production process (e.g. Piran 1999, and references therein), vying against models 
invoking relativistic flows and external shocks for which there is strong Lorentz contraction of the 
timescales between the production mechanism and the observer (e.g. Dermer, Bottcher & Chiang 
1999). Fenimore & Ruiz (1999) have recently argued that a central machine is favored (see also 
Heinz & Begelman 1999). 

One of the implications of this uncertainty is the location of the 7-ray burst. If an external 
blast wave is involved, then a 7-ray burst with a time scale in the observer frame of t Q b s has a 
propagation distance of r 2 ct b s . For a 7-ray burst of 10 s this distance is about 3 x 10 15 cm for 
a Lorenz factor of T ~ 100. This is much larger than the radius of the hydrogen-deficient stellar 
core we are considering here, R CO re — 10 10 cm. It is not clear why the energy produced by shocks 
and UMHDW/LAEMW would be dumped at a radius as large as 10 15 cm. Shocks should deposit 
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their energy as they emerge from the star, and the density falls below the Goldreich-Julian density 
at only 10 11 cm even for a relatively dense wind. If relativistic protons are generated by shocks or 
other mechanisms, they could also plausibly be stopped near the stellar surface by a substantial 
wind. Thus for all the mechanisms we sketch here, the most logical site of the production of any 
7-ray burst is near the stellar surface. This implies that if this general process of pulsar spindown 
from massive star core collapse has any role in producing 7-ray bursts it will most plausibly serve 
as a "central engine." We note that in this case, the natural time scale for any 7-ray burst is about 
5 to 10 sec, the cooling, spin down time for the neutron star. On the other hand, the eruption of 
shocks and UMHDW/LAEMW from the stellar surface may occur in a region small compared to 
the stellar surface, so considerably shorter timescales might be manifest for central burst peaks. 
If the emission from the pulsar is prolonged, then one might also witness time scales associated 
with the processes of production of the UMHDW, for instance various instabilities, as well. The 
shortest time scales of any substructure would be ~ 10 11 cm/T 2 c or about 0.3 ms for V = 100. 

The deposition of a large amount of energy at the stellar surface could, of course, result in 
a subsequent relativistic blast wave and associated afterglow. The processes we are discussing 
would produce a maximum "isotropic equivalent" energy of AitE/AQ ~ 10 54 erg for £52 — 1 
and Ar2/47r ~ 0.01. The external mass required to decelerate this energy, ET~ 2 c~ 2 , requires a 
spherically-distributed mass of ~ 5 x 10~ 5 £52 M Q for T = 100. The mass in the wind out to a 
radius R is 

M wind = 3 x 10~ 6 M & M^v^R 15 , (62) 

with R in units of 10 15 cm. Thus energy emitted by the processes we have outlined will not be 
decelerated until a radius of about 10 16 cm. Even ignoring issues of entrainment and asymmetries 
in the wind density profile, this radius is sensitive to the energy of the burst, the degree of 
collimation, the value of the Lorentz factor and the parameters of the wind. 

The point of this paper is not to establish that core collapse and pulsar formation will 
lead to 7-ray bursts, but to establish that this environment gives a framework in which to 
quantitatively address questions of physics that are germane to the nature of the core collapse 
process and to potential 7-ray production. It seems very clear that rotation and magnetic fields 
have a strong potential to create axial matter-dominated jets that will drive strongly asymmetric 
explosions for which there is already ample observational evidence in Type II and Type Ib/c 
supernovae, their remants, and in the pulsar velocity distribution. The potential to also create 
strong flows of UMHDW/LAEMW serves to reinforce the possibility to generate asymmetric 
explosions. These asymmetries will affect nucleosynthesis and issues such as fall-back that 
determine the final outcome to leave behind neutron stars or black holes. In addition, the presence 
of matter-dominated and radiation-dominated jets might lead to bursts of 7-rays of various 
strengths. The issue of the nature of the birth of a "magnetar" in a supernova explosion is of great 
interest independent of any connection to 7-ray bursts. Highly magnetized neutron stars might 
represent one out of ten pulsar births. Production of a strong 7-ray burst might be even more 
rare. 
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